A formal framework for a nonlocal generalization of Einstein's theory of gravitation 
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The analogy between electrodynamics and the translational gauge theory of gravity is employed 
in this paper to develop an ansatz for a nonlocal generalization of Einstein's theory of gravitation. 
Working in the linear approximation, we show that the resulting nonlocal theory is equivalent to 
general relativity with "dark matter". The nature of the predicted "dark matter", which is the 
manifestation of the nonlocal character of gravity in our model, is briefly discussed. It is demon- 
strated that this approach can provide a basis for the Tohline-Kuhn treatment of the astrophysical 
evidence for dark matter. 
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I. INTRODUCTION 

In special relativity theory, the principle of locality has 
been the subject of a detailed critical analysis and nonlo- 
cal correction terms have been proposed that are induced 
by sufficiently high accelerations This nonlocal spe- 
cial relativity has been applied to electrodynamics and 
the Dirac equation. For a discussion ofpossible experi- 
mental tests of the theory, we refer to [l| and the refer- 
ences cited therein. 

The next subject to be addressed from this nonlocal 
point of view is the theory of gravitation. Einstein's the- 
ory of general relativity (GR) is, by means of the equiv- 
alence principle, heuristically deduced from special rela- 
tivity. Therefore, if special relativity is generalized to a 
nonlocal theory, then general relativity theory cannot be 
exempt from this generalization process. But how should 
such a generalization be implemented in Einstein's the- 
ory? Should one, for example, require the connection to 
be a nonlocal expression in terms of the metric? No ob- 
vious method suggests itself for a direct nonlocal gener- 
alization of Einstein's theory. The underlying reason for 
this is that the precise mathematical form of Einstein's 
principle of equivalence in GR is strictly local; therefore, 
it cannot be employed together with nonlocal special rel- 
ativity to arrive at nonlocal gravitation. On the other 
hand, various heuristic arguments have been advanced 
in favor of a nonlocal classical theory of gravitation (see, 
for instance, @]). 

The idea behind this paper exploits the analogy be- 
tween gauge theories. On the one hand we know that 
electrodynamics, in the framework of a gauged Dirac sys- 
tem, can be understood as a gauge theory of the U(l) 



group; see O'Raifeartaigh 0]. On the other hand it is 
known that a gauge theory of the translation group, for 
spinlcss matter, yields a tclcparallclism theory of gravity 
that, for a suitably chosen Lagrangian, is equivalent to 
Einstein's theory; see, for instance, Nitsch et al. [!, 0]. 
Consequently, if we cannot recognize a direct method of 
generalizing Einstein's theory, it may be helpful to start 
from this so-called telcparallcl equivalent of general rela- 
tivity (GRj|) instead. This will indeed provide a way to 
proceed to a nonlocal extension. 

The analogy between these groups is pointed out in the 
rest of this section and the resulting ansatz for nonlocal 
gravitation is explained in Sec. II. To simplify matters, 
we work in the linear approximation and work out the 
linearized nonlocal gravitational field equations in Sec. 
III. In Sec. IV some subtleties are discussed in connec- 
tion with the use of variational principles in nonlocal 
field theories. Section V is devoted to an interpreta- 
tion of these nonlocal equations in terms of the stan- 
dard linearized general relativity theory but in the pres- 
ence of "dark matter". The properties of the resulting 
"dark matter" are briefly pointed out. A discussion of 
our results is contained in Sec. VI. Various mathematical 
details are relegated to the appendices. In this paper, 
spacetime indices run from to 3 and units are chosen 
such that c = 1. The Minkowski metric tensor is given 
by diag(l, —1, —1, —1). As indices, Latin letters indicate 
holonomic coordinate indices, while Greek letters indi- 
cate anholonomic Lorentz-frame (tetrad) indices. 



A. Electrodynamics 

The electromagnetic excitation Tt^ = (D,H) = —H? 1 , 
a contravariant tensor density, and the field strength 
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with J % as electric current density. Here i, j, ... — 0, 1, 2, 3 
are (holonomic) coordinate indices. The homogeneous 
equations ([IJh can be solved by the potential ansatz 



F« = 2duA n 



(2) 



The two Maxwell equations have to be supplemented by 
a constitutive law that relates the excitation to the field 
strength. In conventional vacuum electrodynamics, one 
takes the local, linear, and isotropic 'constitutive law' 



F 



kl 



(3) 



This law is equally valid for the flat Minkowski spacetime 
of special relativity as well as for the curved Ricmannian 
spacetime of general relativity. In the former case, in 
Cartesian coordinates, the metric gij is constant, in the 
latter case it becomes a field governed by the Einstein 
equations. If we substitute Eq. ([3]) into Eq. |T])i, then 
the Maxwell equations are expressed in terms of the field 
strength i 7 ^ alone. 

In a nonlocal theory, on the other hand, equation (J3]) 
is generalized such that the nonlocal constitutive law ex- 
presses H 1 ^ as a nonlocal spacetime relation in terms of 
Fki- In fact, it is possible to give an effective interpreta- 
tion of Mashhoon's nonlocal electrodynamics in this way 
[f|. That is, one refers the excitation H al3 = e^e^Ti 11 , 
the field strength F a p = e 1 a e 3 pFij , and the metric 
9af3 — e % a€?/3 9ij to non-inertial frames e a = e l a di that 
characterize the tetrads of accelerated observers. With 
respect to such frames, the Maxwell equations read @, 0| 



fyaFp*,] 
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Here C a fp := 2e % a e> pd^e^ is the object of an- 
holonomity ("structure constants") measuring the non- 
inertial nature of the frame. Thus in this approach to 
nonlocal electrodynamics, the nonlocal expression that 
relates the excitation to the field strength is given by @ 

n af3 (x) = V=99 ay 9 PS F lS (x) 

+ [ X a ^ S (x,x')F yS (x')d 4 x'. (6) 



Eq. (|6|) represents a new axiom that, together with the 
Maxwell equations ([1]), determines a nonlocal theory of 
electrodynamics. 



B. Translational gauge theory of gravity 

Let us now turn to the gauge theory of the transla- 
tion group, which includes GR||. We take the field equa- 
tions from || Q ; however ■ we change some conventions 
in order to conform with 0, [l(| ■ For the tensor calculus 
background one should compare Schouten (Tlj . In Ap- 
pendix A we provide a brief sketch of translational gauge 



Expression for 


Electrodynamics 


Transl. Gauge Theory 


conservation law 


di T = 


diTj ~ 


inn. field eqs. 


djTL 3 = J 


djW 3 a ~ T a l 


force density 


fi — F%j ^J^ 


fi ~ Cij a To 3 


horn, field eqs. 


9[iF jk] = 


d[iCj k ] a = 


potentials 


Fa = 2d {l A 3] 


dj a = 2d[iejf 



TABLE I: Analogy between electrodynamics and gravitody- 
namics. The approximation sign employed in the last column 
indicates validity only in the linear regime. 



theory formulated in terms of the calculus of exterior dif- 
ferential forms. Recent works on teleparallelism can be 
found, amon gst others, in the articles of Aldrovandi et 
al. [lj], Itin [I3J3LMaluf et al. [HJll, Mielke [13], 
Obukhov et al. [11 M HI, So et al. [Sf, and Tunget 
al. (22| ; moreover, the two monographs of Blagoievic |23l| 
and Ortin [24[ should be consulted as well as the refer- 
ences cited therein. The inhomogeneous and the homoge- 
neous gravitational field equations are respectively given 
by 



;) yi _ f ' ' 7 ' 



duC. 



0. 



(7) 



The star over the equal sign means that the correspond- 
ing equation is only valid in a suitable frame; this will 
be discussed in Sec. II. The nonlinear correction terms in 
Eq. (JT]) 1 represent the energy- momentum tensor density 
of the gravitational gauge field, 
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As source, we have on the right-hand side of the inho- 
mogeneous field equation the energy-momentum tensor 
density of matter. In analogy with Eq. ((2]), the field 
strength Cij a is defined in terms of the potential e^" by 



2d u 



(9) 



so that the homogeneous gravitational field equations 
02 are thus identically satisfied. 

Clearly, there is a similarity between the electrody- 
namic and gravitational cases, as indicated in Table I. 
In electrodynamics we have one potential Ai as covec- 
tor, corresponding to the one-parameter U{1) group; in 
the gravitational case, because we have four linearly in- 
dependent translations in spacetime, we have four covec- 
tors as potentials e,", the tetrad components with a = 
0,1,2,3. Correspondingly, whereas we have six compo- 
nents of the electromagnetic field strength Fij — —Fij , we 
have 4x6 components of the gravitational field strength 
C i:j a = -Cji a . If one compares Eq. Q with Eq. ©, 
one recognizes the quasi-Maxwellian nature of the grav- 
itational case, only that we have four times more field 
equations in Eq. (Jjj and that the gravitational energy- 
momentum tensor density £ a l emerges because all phys- 
ical fields including gravity are gravitationally 'charged', 
that is, all fields carry energy-momentum. 
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The analog of Eq. ([3]), the local, linear, and isotropic 
constitutive law in GR| | is 0, Q 



H ij a = -( \& 3 a - C a M + 2e [l a C 3] ' 



(10) 



where e := det(ei"). Here k = 8irG, where G is Newton's 
constant of gravitation. If Eq. (flT))) is substituted into the 
inhomogeneous field equations (J7])i, the resulting equa- 
tions have been shown to be equivalent to the Einstein 
equations Taking the electromagnetic case © as a 
prototype, we tentatively expect that a suitable nonlo- 
cal generalization of the constitutive relation (fit))) would 
lead to a nonlocal generalization of Einstein's theory of 
gravitation. It is interesting to contemplate the nature of 
this analogy, which is the basis of the present paper. The 
linear and possibly nonlocal clcctrodynamic constitutive 
law is generally valid for sufficiently weak electromag- 
netic fields, since the constitutive tensor is assumed to 
be independent of the field strength. Therefore, we ex- 
pect that the same holds in the linear approximation for 
the nonlocal theory of gravitation that we develop below 
on the basis of this analogy. 

Finally, we must mention that this approach to nonlo- 
cal gravitation is essentially different from other nonlocal 
modifications of general relativity; see, for example, [25| . 
1201. and the references cited therein. 



II. ANSATZ FOR A NONLOCAL THEORY OF 
GRAVITY 

The arena for teleparallclism is the Weitzenbock space- 
time in which we have the orthonormal tetrad (coframc) 
-d a = ei a dx l . Then the local (anholonomic) metric is 
Va/3 = diag(l, — 1, — 1, — 1) and we can determine the 
(holonomic) coordinate components of the metric by 



9ij — ^lape-i &j 



(ii) 



Simple algebra yields e := det(ei Q ) = ^p^g. 

If a coframe is a coordinate frame, then the object of 
anholonomity 



C a := M° 



(12) 



vanishes; in general, C a is a 2-form which decomposes 
according to C a = \Cij a dx l A dx j such that Eq. (fT2")l in 
components results in Eq. ©. 

The connection 1-form T al3 = -Y^ a = T l al3 dx % in a 
Weitzenbock spacetime is teleparallel. That is, we can 
introduce a suitable global Cartesian tetrad frame with 
respect to which the connection vanishes 



pa/3 JL q 



(13) 



Throughout the rest of our paper we will work in a global 
frame that obeys Eq. (TT3"|) . Then all covariant deriva- 
tives reduce to partial derivatives and we can simplify 



our work considerably. In particular, the torsion of the 
Weitzenbock spacetime becomes 

T a ._ D{ja = M a + Tf3 a ^ ^ ^ = pa ^ ^ 

or, equivalcntly, 



Ti a 



k „ a 



2D [i e 3 f = 2d [l e 3] a -2T [l] fe k 



— 2<9[jej] — dj 



(15) 



Accordingly, the gravitational field strength, in the spe- 
cial "gauge" (TIB"]) , is represented by the object of an- 
holonomity C a . This concludes the geometrical set-up 
of a Weitzenbock spacetime and we will henceforth drop 
the star over the equal sign. 

We now turn to physics in GR| | , which is determined 
via the variation of the action SS = 0, where 



S = J {C g + C m )d 4 x. 



(16) 



Here C g and C m are respectively the gravitational and 
matter Lagrangian densities. Out of the field strength 
C a , we can construct a gravitational Lagrangian £ g that 
is — in analogy with electrodynamics — quadratic in the 
field strength. The explicit form of £ g can be left open for 
the moment. However, we can introduce quite generally 
the components of the excitation Ti 13 a that are related 
to those of the field strength by 



n i3 \ 



(17) 



1 Hij a dx l A 

^e ijM H k i a , where e ljkl is the totally 



In terms of differential forms we have H a — -^±± tja 



dx 3 , with T-L l3 a 
antisymmetric Levi-Civita symbol with e 0123 = 1. Be- 
cause of the presumed quadratic nature of L g and the 
Euler theorem on homogeneous functions, we find 



1 



dG„ 



1 



7~t J uGi 



(18) 



which expresses the gravitational Lagrangian in terms of 
the excitation Ti 13 a and the field strength Cif 1 such that 
TL 13 a is linear in Cij a . 

The field equations of GR | given in Eq. ([7]) follow from 
the variational principle of stationary action (|16p with 
T a l := 5C m /Sei a as the source. We note that, in lin- 
ear approximation, the gravitational energy-momentum 
complex ([5]) will vanish and only the first term of the 
left-hand side of Eq. ©i will survive. 

Having thus established the geometry and the 
Lagrange-Noether machinery for the teleparallel theory 
of gravitation, we turn next to the explicit choice of the 
gravitational Lagrangian. General relativity is recovered 
via equation (| 10|) . which may be written as 



(19) 
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with the modified field strength 
1 



(20) 



In particular, along the same line of thought as in 
Q, we wish to consider a nonlocal constitutive law for 
gravitation. To prepare the ground, let us introduce the 
invariant proper infinitesimal distance ds between two 
neighboring events in Wcitzcnbock spacetimc 



lis 2 = gijdx 1 (g> dx^ 



(21) 



and define a geodesic between two fixed events P' and 
P to be the path that is an extremum of the spacetime 
distance between P' and P, 



5 J ds = Q. 



(22) 



This path is given by the geodesic equations 

dPx 1 ( i } dx j dx k „ , . 

where 



jkl = \ ( " 9li ' k + mk ' 3 ~ 9jk ' 1 ^ ^ 



are the Christoffel symbols. 

We assume that two causally separated events are con- 
nected by a unique timelike or null geodesic; more gener- 
ally, in the spacetime region under consideration, there 
exists a unique geodesic joining every pair of events. 
It then proves useful to employ the world-function O, 
which denotes half the square of the proper distance from 

P' : x' = £(Co) to P : x = £(Ci) along the geodesic path 
x i = That igj we define m 

1 dP dP 

0(x,xO = -(Ci-Co)^ dij^-^dC (25) 

It turns out that O is independent of the affinc parameter 
C; moreover, the integrand in Eq. (|25[) is constant by 
Eq. (|23j) . The main properties of Cl(x, x') are summarized 
in Appendix [Bl 

To distinguish coordinate indices that refer to x from 
those that refer to x' , we henceforth use indices a, b, c, ... 
at x and k, ... at x' . Thus we define 

,90 , dtt 

Q ^ x )=g x ^> Q i( x > x ) = 7)^i> ( 26 ) 



dx' 



and note that covariant derivatives at x and x' com- 
mute for any bitensor. It follows from the results of Ap- 
pendix [B] that 



Differentiating this equation, we find that Q a i(x, x') = 
VL ia (x,x') are smooth dimensionless two-point tensors 
such that 



lim to ai (x, x') = -g a i{x) . 

X 1 — >X 



(28) 



Thus a possible nonlocal generalization of Eq. (fl9|) is 
given by 



n ab c (x) 



V^hlxj / u(x,x')n ai n b m ckX (x,x') 



x£ ij k (x') y /-g(x>)d i x' , 



(29) 



where U(x, x') is unity except when x' is in the future 
of a;, in which case U vanishes; in Minkowski spacetime, 
this means that events with a;' > x° are excluded from 
the domain of integration in Eq. (|29p . Moreover, x( x , %') 
is a scalar given, for instance, by 

X{x,x') = 5(x-x') + K(x,x'), (30) 

where the Dirac delta-function S(x — x') is defined via 



5(x - x')ip(x') ^f-g{x') d 4 x' = ip(x) (31) 



for any smooth function ip(x). The scalar kernel K(x, x') 
in Eq. (|30|) denotes the nonlocal deviation from a local 
constitutive law. The main nonlocal constitutive relation 
([29"]) can therefore be expressed as 

1 

K 



n ab c (x) = -^gjx)[€ ab c (x)- J u(x,x')a ai n bj n ck 

xK(x, x , )£ ij k (x')y/-g(x') d 4 x'} , (32) 



where the nonlocal deviation from general relativity is 
made explicit. More complicated nonlocal constitutive 
relations are certainly possible; however, Eq. (|32[) is the 
simplest one involving an unknown scalar kernel K. 

We emphasize that the domain of applicability of 
Eq. (|32p is not physically restricted; in particular, a 
nonlocal gravity theory of this type is not directly re- 
lated to nonlocal special relativity [l[. That is, if the 
whole heuristic machinery of Einstein's principle of equiv- 
alence could be employed in this case, one would have 
to conclude — on the basis of nonlocal special relativity — 
that gravity should be nonlocal for sufficiently high gravi- 
tational "accelerations". However, Einstein's principle of 
equivalence cannot be used in this way, as its precise for- 
mulation in GR is the very embodiment of locality. Thus 
nonlocal gravity is in this sense decoupled from nonlocal 
special relativity. 

It is important to note that our ansatz (|32p is non- 
linear as well as nonlocal, since, among other things, the 
metric tensor (|11[) is quadratic in the gravitational poten- 
tials ei a . The scalar kernel K(x,x') could be given, for 
instance, by K~ x = 2 + Lq, where Lq is a constant char- 
acteristic length. It could also involve the Weitzenbock 
invariants Ial 



20 = g ab n a n b = gVfliClj 



(27) 



^ijk^ i ^kji^ j ^ij ^ k 



(33) 
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at x and x' as well as scalars formed from the covariant 
derivatives of f2(x,x'). To illustrate the latter, consider, 
for example, fl a and fij in Eq. ([26]) . The vector fi a is 
tangent at x to the geodesic connecting x' to x and the 
length of il a is equal to that of this geodesic by Eq. (f27]k 
in Minkowski spacetime, £l a = x a — x la . Similarly, the 
vector f2 z has the same length, is tangent to the geodesic 
at x', and is directed away from x; in Minkowski space- 
time, ft 1 = x n — x l . Hence one can form coordinate 
scalars £l a e a a {x) and fi*ej a (x'), which turn out to be of 
particular significance in Sec. V. This brief discussion in- 
dicates that many options indeed exist for generating a 
scalar kernel K(x, x') in our ansatz ([32]) . 

In the following section, the nonlocal theory of grav- 
ity is treated in the linear approximation; therefore, 
Eq. (|32p will be employed in linearized form. In par- 
ticular, U(x,x')K(x,x') will be evaluated in Minkowski 
spacetime, since Weitzenbock spacetime in the special 
gauge (fT3|) reduces to Minkowski spacetime when the 
gravitational potentials e^" reduce to 5f. In this limit, 
we define 



K (x, x ) = U (x, x')K(x, x 1 ) 



(34) 



III. 



GRAVITATIONAL FIELD EQUATIONS IN 
LINEAR APPROXIMATION 



The gravitational field equations are obtained from 
the variation of the action, J £d 4 x, where C is the La- 
grangian density given by the sum of the correspond- 
ing quantities for gravitation £ g and matter £ m , namely 
C = C g + C m , such that 



Sei a 



(35) 



Here T a l is the energy-momentum tensor of matter, so 
that y^gT a l = T a l given in Eq. (J7J). The gravitational 
part £ g is given by Eq. (fl~8|) in analogy with electrody- 
namics 0- For the nonlocal theory under consideration, 
a is now given by Eq. (|32p ; however, to simplify mat- 
ters, we will work in the linear approximation. That is, 
we assume that 



(36) 



where the nonzero components of ip % a are such that 
\'4 >% a\ <C 1. In this approximation, the holonomic co- 
ordinate indices and the anholonomic tetrad indices are 
indistinguishable. It then follows from Eq. (|32[) and Ap- 
pendix [B] that 



1 

where K is the scalar kernel in the weak-field approxima- 
tion. A discussion of action principles in nonlocal field 



^<V(x) / K(x,y)^ k (y)<fy, (37) 



theories is contained in Edelen [28j. The expressions for 
Cij k and C y fc in terms of can be easily determined. 
We find from Eq. ([H that 



a j:a 



so that from ©, 



S? + il> a i 



Moreover, 



It follows from (gOD that 



(38) 
(39) 
(40) 



-- (h i j - h 3 
2 V k ■ k 

+ 6i(t j -i> l j t l )-5i(t i -Tpi i , 1 ) , (41) 



where i\) = rjij ip v ■ 

Let us first examine the local part of the Lagrangian 
(l37l). It follows from 



<J£ g (local) = (SC t] k ) - i- (8& k ) Cij k (42) 



(43) 



that the first term in Eq. (|42|) is in effect given by 



via Eq. (|39p and integration by parts, assuming that Sip k i 
vanishes on the boundary, and hence neglecting such 
boundary terms in Eq. (|43p . Applying this same pro- 
cedure to the second term in Eq. (|42|) . we find after a 
detailed calculation that the result is again given by 



2 k 



(44) 



so that in the absence of the nonlocal term [K = 0) 
the variation of the action can be expressed in the linear 
approximation as 

( + TV 1 ) (S^i) d 4 x = . (45) 

Thus for K = 0, the gravitational field equations are 
dj<£f 3 k = nTk 1 ; moreover, we prove in Appendix ICl that 
these are precisely Einstein's field equations according 
to the standard general linear approximation scheme. A 
general proof,_yalid in the non-linear regime, is already 
contained in Q. We note that 3j2V = is implied by 
the field equations. 

We now turn to the nonlocal part of the Lagrangian. 
In this case, the corresponding two terms in the variation 
of the gravitational action are 

<5£ g (nonlocal) = ~ [SCij k {x)] j K(x,y) €h(y)d A y 

-~k Cijk{x) I K{Xj y) [^ fc(2/) ] d ' y • (46) 
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As before, the first term in Eq. (|46| is easily shown to 
lead to 

-^[WV*)] J^^^^y. (47) 

For the second term in Eq. (|46[) . we can write its contri- 
bution to the variation of the total action as 

" hi I [8&k(yj\ d 4 yd 4 x. (48) 

The domains of integration in the double integral are the 
same; therefore, it is possible to switch x and y in Eq. 1(15]) 
to get 

~hj j K ^ x ) c ^v) i^kW] d vv ( 49 ) 

Hence the second term in the nonlocal part of the La- 
grangian is equivalent to 

- i- [<^' fc (*)] | (y, x)C ij \y)d i y . (50) 

Now applying to Eq. (|50[> the same procedure we used 
in the derivation of Eq. (H3}, we find after a detailed 
calculation that the result is 

Let us define the symmetric kernel K by 

K(x,y) = ^[K(x,y) + K(y,x)] . (52) 
Then, combining equations (|47j) and (|51|). we find 
<5£ g (nonlocal) = [t5?/' fc i (cc)] 

x/%^ fc (^V (53) 

This means that the field equations of the nonlocal theory 
in the linear approximation are given by 

dj&k + J dK Q X J V) &h{v)#V = «ZV . (54) 

We note that dCT^ — is still identically satisfied in 
the nonlocal theory. For a general symmetric kernel K, 
the energy-momentum tensor Ty is not symmetric in gen- 
eral, since the nonlocal part in Eq. (|54[) is not in general 
symmetric. This poses no basic difficulty as Ty, given 
in general by Eq. (|35p . is not symmetric in general. On 
the other hand, the requirement that IJy be symmetric 
(as in GR, for instance) would impose restrictions on the 
kernel K. 

A natural interpretation of the basic nonlocal gravita- 
tional field equations (|5^|) can be obtained by consider- 
ing the nonlocal term to be an effective source for "dark 



matter" and hence moving it to the right-hand side of 
Eq. (|5^|) . In this way, one can recover Einstein's theory 
but with "dark matter". We explore this possibility in 
the rest of this paper. 

Consider, for instance, the possibility that K is an even 
function of x — y; that is, 

K(x,y)=F(x-y), (55) 

where F(z) = F(—z). Then, dK/dx^ = —dK/dyi and 
Eq. ([Ml) can be written as 

dj& k (x) + J K(x, y) ^M d 4 y = KTj , (56) 

which implies, via Appendix C, that Tki = Tik- Here we 
have assumed that the derivatives of ip k i vanish on the 
boundary of the spacetime region under consideration. 



IV. NONLOCAL FIELD EQUATIONS 

The treatment of the previous section has been based 
on the assumption that nonlocal gravitational field equa- 
tions should be obtained from an action principle that 
incorporates the nonlocal constitutive relation (|32|) . Al- 
ternatively, this relation could be employed in the field 
equations of tclcparallelism. The purpose of this section 
is to address certain subtleties involved in these possibil- 
ities. 



A. The "general" field equations of teleparallelism 

In electrodynamics the Maxwell equations in fTj can 
be derived from electric charge and magnetic flux con- 
servation, see Q. No action principle is necessary. If 
one considers, for example, local and linear magnctoclcc- 
tric matter, the constitutive law is = oX^Fkl, with 
a constitutive tensor density x l ^ kl = ~x = ~x 
that has 36 independent components. This model can 
be derived from a Lagrangian, provided the constitu- 
tive tensor obeys additionally the Onsager-type relation 
yijki _ ^khj . -j-jjg^ jf £ 3g components vanish. 

Nevertheless, even if this condition is not fulfilled, one 
can still find reasonable applications in physics, though 
in this case irreversible processes have to be taken into 
account. 

The situation is similar for the energy-momentum ten- 
sor density of the electromagnetic field. Starting from the 
Lorentz force density, substitution of the inhomogeneous 
Maxwell equations and partial integration result in the 
(canonical) Minkowski energy-momentum tensor density 

7? = -^> H « W + F ik W k . (57) 

Apparently no Lagrangian is necessary for the derivation 
of this expression 0, [29| . 
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The gravitational field equations of the translational 
gauge theory, as formulated in Eq. 0, should also have 
general validity within this framework, independently 
of the existence of the Lagrangian. For the homoge- 
neous equations 02 this is apparent since they represent 
just the first Bianchi identities in a Weitzenbock space- 
time. For the inhomogencous equations ([7])i this can 
be seen as follows: We start with the energy-momentum 
conservation law for matter in the linear approximation 
d-Tj ~ 0. We can "solve" it by T a l = djW* Q - Ej and 
TC 10 a = —H. Jl a , with some unknown correction term . 
These are already the inhomogeneous field equations. We 
just have to determine E a % explicitly. Clearly it adds to 
T a l , namely djH l] a — T a l + E a % \ then, it is natural to 

interpret E a l as the energy-momentum tensor density of 
the gravitational field. In analogy with Eq. (|57p we ex- 
pect it to have the form of Eq. ([5]) ; we now drop the tildes 
from E a % and H 10 a . Therefore the inhomogeneous field 
equations extracted from the energy-momentum conser- 
vation law of matter in this heuristic manner have exactly 
the same form as (JT]) 1 with |(5J). Thus we may postulate 
([7]) as the general field equations valid independently of 
the existence of a Lagrangian. 

By contrast, if we have a Lagrangian, the field equa- 
tions turn out to be § 

djH ij a - (e\C s + Cj'H^fs) = T a % . (58) 

If £ g = — jCij a H l: > a is substituted into l[5"8"|). we recover 
0i- Consequently, we have shown that the "general" 
field equations 0i with ([5]) are correct if a Lagrangian 
exists. Accordingly, they are certainly one consistent and 
reasonable generalization of the Lagrangian-bascd equa- 
tions (|55|) . 

B. Ambiguity in the field equations 

We now wish to address a certain ambiguity that is 
encountered in implementing a nonlocal constitutive law 
as in our present work. A local constitutive law relating 
the excitation "W 3 k to the gravitational field stregth C %3 \ 
can be used in the derivation of the field equations of the 
theory in cither of two equivalent ways: it could be em- 
ployed in the Lagrangian (TT5|) that then results, via the 
variational principle of stationary action, in the desired 
field equations or, alternatively, it could be directly sub- 
stituted in the field equations 01. These two methods 
in general produce different results for a nonlocal con- 
stitutive law, however. It is interesting to illustrate this 
point for the case at hand in the linear approximation. 
In this regime, we have the linear constitutive law 

KH l] a {x) = <£ ij a (x) + J K(x,y)C j a (y)d 4 y , (59) 

which, when inserted in Eq. (fT5|) results, as described 
in detail in Sec. Ill, in the nonlocal field equations ([51)) . 



Alternatively, however, we can equally well substitute the 
constitutive relation (|59[) in equation 0i, which reduces 
in the linear regime to the new nonlocal field equations 

dj&k + J 9K g% V) &k(v)d*y = KTk* ■ (60) 

The two nonlocal field equations only differ in their ker- 
nels: Eq. (f5"l)) involves K, while Eq. ([60|) involves K. 
The question of whether such nonlocal equations as (|60|) 
can be derived from a variational principle is beyond the 
scope of our investigation; for a related study in connec- 
tion with acceleration- induced nonlocality see [30j ] - 

In a similar way as in Sec. Ill, one can recover a direct 
nonlocal generalization of Einstein's theory with a sym- 
metric energy-momentum tensor of matter by assuming 
that K (x, y) is a function of x — y. In this case, Eq. (|6"0)) 
can be written as 

d 3 €h(x) + J K(x, y) d ^ v) #y = KTj , (61) 

which should be compared and contrasted with Eq. (|56[) . 
These inequivalent nonlocal field equations both admit 
a natural interpretation in terms of dark matter as 
described in detail in the following section; however, 
Eq. (|6"Tj) has an advantage over Eq. (f5"6"]) in that it in- 
volves a causal kernel K, while in the symmetric kernel 
K of Eq. (f56"|) past and future are treated in the same 
manner. 



V. "DARK MATTER" 

Let us now consider an important consequence of our 
nonlocal equations for the gravitational field. To keep 
our discussion general, we start with an equation that 
has the same form as Eqs. (|5^|) and (fBI)]) , except for a 
kernel /C, which we take to be equal to K or K, respec- 
tively, depending on whether one adopts the Lagrangian- 
bascd approach of Sec. Ill or the more direct approach 
of Sec. IV. We henceforth assume, for simplicity, that in 
the former case K(x, y) is an even function of x — y and 
in the latter case K(x,y) is simply a function of x — y. 
Then, the arguments of the previous sections imply that 
the nonlocal modification of Einstein's gravitational field 
equations for a symmetric energy-momentum tensor T^- 
may be expressed as 

Gij(x) + J !C(x,y)G l3 (y)d 4 y = kT^x) , (62) 

where Gij is given by 

G { j = d k € k , (63) 

and represents Einstein's tensor in the linear approxima- 
tion. 

The nonlocal term in equation (|62[) can be interpreted 
in terms of an effective energy-momentum tensor for 
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"dark matter" by writing Eq. ([62)1 as 

Gij(x) = k [Tij(x) + %ij(x)] , 



(64) 



where the "dark matter" component is given by the sym- 
metric energy-momentum tensor 



%ij(x) = -~ / K(x,y)G ij (y)d 4 'y. 



(65) 



Equation (|65|) may be written in a more transpar- 
ent form by making use of the methods of Appendix |D] 
That is, suppose that the integral equation (|6"2")l can be 
solved by the method of successive substitutions (cf. Ap- 
pendix [D]) and that the resulting infinite series is uni- 
formly convergent. Then, as shown in Appendix D, it 
is possible to introduce a reciprocal scalar kernel 1Z such 
that 

Gij(x) = nTij(x) +K Tl(x,y)T ij (y)d i y . (66) 



Thus the "dark matter" component in equation (|64j) is 
given by 

% 3 (x) = J K(x,y)T l3 (y)d 4 y, (67) 

which is the integral transform of matter T^- by the ker- 
nel TZ(x,y). It is clear from Eq. (|6T[) that in our model, 
"dark matter" should be quite similar in its character- 
istics to actual matter; for instance, for dust, the corre- 
sponding "dark matter" is pressure-free, while for radia- 
tion with Tk k = 0, the "dark" energy-momentum tensor 
%ij is traceless as well. 

The reciprocal kernel 1Z is given as an infinite series in 
terms of iterated kernels constructed from IC(x,y); that 
is, 



fc(x,y) = ^2^n{x,y) 



(68) 



To ensure causality, it is useful to assume further that 

K(x,y)=6(x - V )P(x,y), (69) 

where P is a function of x — y. Then, all iterated kernels 
as well as the reciprocal kernel are also of this general 
form; indeed, 



Kn(x,y) = 6(x°-y°)P n {x,y), 
where Pi(x, y) = P(x, y) and 

P„+i(x,y) = -J P(x,z)P n (z,y)d 3 



Therefore, 



K(x,y) = 6(x»-y»)Q(yL,y), 



(70) 



(71) 



(72) 



where Q is reciprocal to P, 



<2(x,y) 



P«(x,y) 



(73) 



Moreover, P n and Q are functions of x — y, since the 
integration in Eq. (|7ip extends over the entire three- 
dimensional Euclidean space. Thus "dark matter" is 
the convolution of matter and the reciprocal kernel in 
this case. The mathematical implications of this fact are 
treated in Appendix E. 

It is interesting to consider the confrontation of our 
nonlocal theory in the linear approximation, represented 
by equations ((62|) - ([67|) . with experimental data. Let us 
first determine the Newtonian limit of our nonlocal grav- 
ity. This follows directly from Eq. (|64[) in the standard 
manner and we find the Poisson equation 



V 2 $ = 47rG(p + p D ) , 



(74) 



where $ is the Newtonian potential. Here p is the density 
of matter and /?d, 



p D (t,x)= J Q(x,y)p{t,y)d 3 y , 



(75) 



is the corresponding density of "dark matter" . We note 
that in the general linear approximation of our theory the 
kernel 1C{x,y) and its reciprocal 1Z{x,y) are given func- 
tions in Minkowski spacetime and are thus independent 
of any particular physical system. The solar-system tests 
of general relativity imply that our "dark matter" must 
be a rather small fraction of the actual matter in the so- 
lar system. This can be simply arranged with a suitable 
choice of the universal reciprocal kernel Q. 

It is interesting to illustrate these considerations in the 
case of the problem of dark matter in spiral galaxies (3ll . 

mm in us mull- 

Imagine, for instance, the circular 
motion of stars in the disk of a spiral galaxy. Beyond the 
galactic bulge, the Newtonian acceleration of gravity for 
each star at radius |x| is given by Vo 2 /l x l toward the center 
of the galaxy. Here is the constant "asymptotic" speed 
of stars in accordance with the observed rotation curves 
of spiral galaxies. It follows from Poisson's equation that 
the effective density of dark matter is essentially given by 
u 2 /(47rG|x| 2 ). Comparing this result with equation ([75]) 
and setting 



p(t,y) = MS(y), PD(t,x) 



4ttG|x| 



(76) 



where M is the effective mass of the galaxy and the di- 
mensions of the galactic bulge have been ignored, we find 



Q(x,0) 



47rGi\/|x| 2 



(77) 



The reciprocal kernel Q is a function of x — y; therefore, 
it follows from Eq. (J77)) that 



Q(*,y) = A 



47rA |x - y| 2 



(78) 
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with a universal length parameter A = GM/vq. Thus 
taking due account of the observed rotation curves of 
spiral galaxies, Eqs. (f74"|) and (|75|) with kernel (|78|) imply 



= 4ttG 



p(*,x) + 



1 

4^A 



p(t,y)d 3 y 
|x-y| 2 



(79) 



It is important to point out in passing a defect of the spe- 
cific form of Eq. (f78|) : the total amount of "dark matter" 
associated with a nonzero matter density p is infinite (see 
Appendix E); therefore, kernel ([78]) is too simple to be 
quite adequate for the task at hand. 

For a point mass m with p(i,x) — mS(x), the Newto- 
nian potential given by Eq. (|79| can be expressed as 



, Gm Gto , / Ixl 



(80) 



The observational data for spiral galaxies indicate that 
A is of the order of a kpc, so that the logarithmic term 
in Eq. (|80|) is essentially negligible in the solar system. 
Moreover, the universality of the kernel in the linear ap- 
proximation implies, via A = GM/vq, that for spiral 
galaxies M on Vq, since A should be independent of any 
particular physical system. 

It is remarkable that in the simple considerations re- 
garding Eqs. (|79|) and (|80|) — based on our linear approx- 
imation scheme — we have recovered the significant pro- 
posal put forward by Tohline and further developed by 
Kuhn et al. [H, [39| to solve the dark matter problem 
by a natural modification of the Newtonian law of grav- 
itation. An interesting discussion of the Tohline-Kuhn 
scheme is contained in the review paper of Bckenstcin 
[39l ] . Despite various successes, the main drawback of 
this approach appears to be the violation of the Tully- 
Fisher relation, which implies that M cx Vq [H, [39[ . To 
agree with the empirical Tully-Fisher law, it seems nec- 
essary to go beyond the linear approximation scheme. 

Within the general framework of this work, but going 
beyond the linear approximation as well as the specific 
constitutive model employed thus far, it may be possi- 
ble in principle to have a relation of the general form 
of Eq. (fr?T)) with a kernel that is highly dependent upon 
the particular system under consideration. This could 
provide a natural way to interpret the observational ev- 
idence for dark matter as the nonlocal manifestation of 
the gravitational interaction. It remains to elucidate the 
physical origin of the constitutive kernel that has been 
the starting point of our investigation. 



VI. CONCLUSION 

To develop a nonlocal generalization of GR, it proves 
useful to approach Einstein's theory via its equivalent 
within teleparallclism gravity, namely, GR||. Therefore, 
we work in Weitzenbock spacctimc with a tetrad field 



ei a and its dual e^p, a metric gij, and a flat connec- 
tion Ti a P = —Ti^ a that is chosen to vanish globally 
(Ti a/3 = 0). The gravitational field strength is given by 
Cij a and the modified field strength by <£jj a . In this 
framework, which is capable of nonlocal generalization, 
GR|| corresponds to a specific gravitational Lagrangian. 
Working with a nonlocal "constitutive" kernel K (x, x') in 
the linear approximation, we construct an explicit nonlo- 
cal generalization of Einstein's theory of gravitation that 
is consistent with causality. This theory can be refor- 
mulated as linearized general relativity but with "dark 
matter" , which mimics the contribution of nonlocal grav- 
ity. We find that the effective energy-momentum tensor 
of "dark matter" is simply the integral transform of the 
energy-momentum tensor of matter. 

The application of our nonlocal model in the linear ap- 
proximation to the dark-matter problem in spiral galax- 
ies is in conflict with the empirical Tully-Fisher relation. 
It is possible that the situation can be significantly im- 
proved with a suitable choice for the kernel in the non- 
linear regime. However, a more basic theory is needed 
to determine the nonlocal constitutive kernel from first 
principles. This is a task for the future. 

Acknowledgments. We are grateful to Jacob Bekcn- 
stein (Jerusalem), Yakov Itin (Jerusalem), Jose Maluf 
(Brasilia), Eckehard Miclke (Mexico City), and Yuri 
Obukhov (Moscow) for helpful comments. 



APPENDIX A: TRANSLATION AL GAUGE 
THEORY IN EXTERIOR CALCULUS 

The coframe 1-form d a = ei a dx l represents the gravi- 
tational potential; here i, j, ... = 0, 1, 2, 3 are (holonomic) 
coordinate indices and a,/?,... = 0,1,2,3 arc (anholo- 
nomic) frame indices; the frame ep = e 3 pdj is dual to 
the coframe, that is, ei a e l p = 5% and ei a e 3 a = S 3 . The 
field strength of gravitation is the object of anholonomity 
2-form C a :=dti a . 

Spacctimc is described by a Weitzenbock geometry 
with a telcparallcl connection 1-form Y a P = Tia^dx 1 and 
with the curvature 2-form 

Rj ■= dTj - r a 7 A r/ = -Rijjdx* A dx 3 (Al) 



that vanishes 



R a p = . 



(A2) 



We decompose the connection into a Riemannian (Levi- 
Civita) part and the contortion, 



r p = f f 3 -k p 



(A3) 



with the following definitions of the torsion 2-form and 
the contortion 1-form, respectively, 

T a := Dd a = dtf a + T p a A d 13 , T a =: K a A d p . (A4) 
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As in a Euclidean space, one can likewise choose in a 
Weitzenbock space suitable frames such that the connec- 
tion vanishes everywhere: 

r</ = TiJda? = . (A5) 

Eq. (|A5j) should be understood as a choice of a certain 
gauge, which is possible since the curvature R a ^ van- 
ishes everywhere. Many calculations are simplified in this 
gauge. 

The Lagrangian 4-form is the sum of a gravitational 
part and a matter part 



L = L s (# a ,C a ) + L m (ip,di>,# a ) , 
The field equation reads 



dH a — E a = 



6$ a 



(A6) 



(A7) 



Here we have the 2-form H a := —dL g /dC a = ^Hij a dx l A 
dx J as the excitation and the 3-form E Q as the energy- 
momentum density of matter, whereas the 3-form E a 
represents the energy-momentum density of the gravita- 
tional field 111 



E a :=e Q jL g +(e Q JC^) A Hp 



(A8) 



Equations (|A7[) with (|A8[) represent the Lagrangian- 
based form of the field equations. If we express the 
gravitational energy-momentum 3-form E a in terms of 
a Minkowski type current [or, equivalcntly, if we substi- 
tute L g = -\C a ATL a into Eq. (|A8[) ], we find the general 
field equations 

dH a - i [Hp A {e a \CP) — C 13 A (e a \Hp)] = S Q . (A9) 

This is the exterior-form version of the field equations 
©i with Eq. ©. 

Usually H a is, like in electrodynamics, linear in the 
field strength 

Ha = M *( fll (DCT/S + a2 (2) C /3 + a3 (3) C /3) . (A10) 

Here ^C' 3 are the different irreducible pieces of C@ and 
* represents the Hodge star. The "Einsteinian" choice 
for ^C 13 , I =1,2, 3, turns out to be 



di = —1 , 0,2 = 2 , 03 



(All) 



In this case, Eq. (|A10|) is distinguished from all linear 
Lagrangians in that it becomes locally Lorcntz covariant. 



On the other hand, 



so that 



dD, . „ dfl . ,, 
dO = — — Sx a + — — -rSx , 

ox a ox" 



on dx» 

= ^-^)9a b (x) — 



(B2) 



an 



dx' 3 



dx a -(Ci-CoK^') rfC 



(B3) 



It is possible to see from Eq. (|2"3"|) that the integrand in 
Eq. (|2"5j) is indeed constant; therefore, 

n(x,x') = i(Ci - MViM^^r 

= ife-O^Mf £ • (B4) 

It follows from Eqs. jB3|-(jB4| that Eq. (J2TJ) is satisfied; 
moreover, f2 = for a null geodesic, = ^ r 2 for a time- 
like geodesic of length r and fi = — i cr 2 for a spacelike 
geodesic of length a. 

In Minkowski spacetimc Q is given by 

Q(x, x 1 ) = ±mj (x H - x l ){x' 3 -x 3 ). (B5) 

According to our convention, r/ij = diag(l, — 1, — 1, — 1); 
hence, 

n(x,x') = \[(t'-t) 2 -(x'-x) 2 ] . (B6) 



In this case, we find that 

d 2 n(x,x') 



dx a dx n 



Vai 1 



while 



d 2 n 

dx a dx b 



Vab , 



d 2 n 



13 ftn-H 



dx H dx' 3 



(B7) 



riij . (B8) 



APPENDIX C: SYMMETRY OF THE TENSOR 



The purpose of this appendix is to show that dj<L^ k is 
a symmetric tensor, so that in the field equations 



dj€ j k = KlV 



(CI) 



APPENDIX B: PROPERTIES OF THE 
WORLD-FUNCTION Q, 

Consider a variation of Eq. ([23)) that changes the end- 
points, then 



«f!(i,i') = (Ci-Co) 



Ci 



0) 



(Bl) 



the energy-momentum tensor of the source is symmetric 
(Tjfc = Tfei). In fact Eq. (|C1[) is identical to Einstein's 
field equations in the linear approximation. 
It follows from Eq. (gj) that 

dj&k = -\{Uh k i -h kj , i3 )+\h i ^ kj 



sih^jt + sin^p-t'k, (C2) 
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where we have used the relation 2ijji\ji — hp ji. Define 
the trace-reversed gravitational potentials hn~ as 



(C3) 



where h = rjijh ZJ = 2ip. Then, replacing hij by hij + r]ijip 
everywhere in Eq. (|C2p results in 



fytjk = -\uhki + \h k \ t] + ±hi' M - \mih Jl 3 , (C4) 

which is manifestly symmetric in i and k. From Eqs. (|C1[) 
and (fC4|) . we get 

-□/life + h/ ijk + h k\ 3 i ~ Vikh J \ji = 2nT lk . (C5) 

These are exactly the same as Einstein's field equations 
in the linear approximation. 



APPENDIX D: LIOUVILLE-NEUMANN 
METHOD OF SUCCESSIVE SUBSTITUTIONS 



Then, Eq. (|D2[) can be written as 



f(x) = <f>(x) + / r(x,y)f(y)dy. (D5) 



It is clear from equations (|Dip and (|D5|) that the kernels 
k and r are reciprocal of each other. 

Let us note here some properties of these kernels. It 
follows from Eq. (|D3|) that 



k n+p (x,y)= / k n (x,z)k p (z,y)dz , (D6) 



where p = 1,2,... . Using Eqs. (|D3j) and (|D4|) we find 
that 

,6 

k(x,y) + r(x,y) = / k(x, z)r(z 7 y)dz 

J a 
fb 

r{x, z)k(z,y)dz . (D7) 

Moreover, if k is symmetric, k(x,y) = k(y,x), then all 
iterated kernels as well as r(x, y) are likewise symmetric. 



Consider a linear integral equation of the second kind 
given by 



x )=f( x )+ k(x,y)(j)(y)dy , 



(Dl) 



where a and b are constants. We seek a solution of this 
Frcdholm equation by the method of successive substitu- 
tions. That is, we replace <f> in the integrand by its value 
given by equation (|D1|) . Repeating this process eventu- 
ally leads to an infinite series of the type 

<t>(x) = f(x)+ / k(x,y)f(y)dy 

J a 

k{x,z)k(z,y)f(y)dydz + -- - . (D2) 

If this series is uniformly convergent, we have a solution 
of the integral equation (|D1|) . This solution is unique 
in the space of real continuous functions on the interval 
[a, 6]; a generalization of this result to the space of square- 
integrable functions is contained in Tricomi [io| . 

Let us define the iterated kernels k n , n — 1,2, ... , by 

h(x,y) = k(x,y), 

k n+1 (x,y) = / k(x,z)k n (z,y)dz . (D3) 



These functions occur in the infinite series of Eq. (|D2 
We define the reciprocal kernel r(x,y) such that 



r(x,y) = ^2,k n {x,y) . 



(D4) 



APPENDIX E: CONVOLUTION KERNELS 

Suppose that the spatial kernels in Sec. V are such that 

P(x,y)=p(x-y), Q(x,y) =«(x-y). (El) 

The purpose of this appendix is to point out the conse- 
quences of the convolution theorem for a system of den- 
sity p such that the density of "dark matter" pd is a 
convolution of p and q. In general, the spatial kernels 
p and q could depend on the characteristics of the par- 
ticular system under consideration. The main results of 
interest here are Eqs. (|69[) to (|75|) . Let us note that the 
particular solution of Poisson's equation (|74[) is given by 



$(t,x) = -G J 



p(t,y) +pu(t,y) 3 

|x-y| 



d?y-, 



(E2) 



provided $ vanishes sufficiently fast at spatial infinity. It 
follows from Eq. (|75|) that this Newtonian potential can 
also be expressed as 



where 



*(t,x)=-G / 5(x-y)p(t,y)d 3 2/ , 



q(z)d 3 2 



5(r) = R + / 



(E3) 



(E4) 



Let us assume that the functions of interest here can 
be expressed as Fourier integrals; that is, for a function 



vf{v) = / f(k)e lk - T d 6 k 



(E5) 
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where /(k) is given by 

i//(k) = / f{r)e~ lk r d 3 r . 



It follows that 



(E6) 



For the sake of simplicity, we have introduced here a con- 
stant parameter v, 



(27T) 3 / 2 . 



(E7) 



We are interested in Fourier integral transforms of real 
functions; therefore, /*(k) = /(— k). 

It follows from the convolution theorem for Fourier in- 
tegrals that Eq. ([75]) implies 



pbQn) = vp(k)q(k). (E8) 
Here q(k) is a dimensionless function such that 



vq{0) = ~M 



(E9) 



where M = ^p(O) is the mass of the system under con- 
sideration and Md = ^pd(O) is the corresponding "dark" 
mass. Consider, for instance, kernel {75} that is associ- 
ated with the discussion of the rotation curves of spiral 
galaxies in Sec. V; that is, 



3(r) 



1 1 

4^A W 



(E10) 



i i 



4(2tt) 1 /2A |k| 



(Ell) 



Thus for any M > 0, M D = oo. 

Once g is determined from Eq. (|E8|) . it is possible to 
work out its reciprocal kernel. In fact, it follows from 
Eqs. {7_D and {73} that p x (k) = p(k), 



and 



p n +i( k ) = ~^P( k )P«(k) , 



(k)=£p n (k). 
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It is then straightforward to show that 

p(k) 



Therefore, 



1 + vp{k) 
1 + uq(k) 



(E12) 



(E13) 



(E14) 



(E15) 



These results are consistent with the fact that p and q 
are reciprocal of each other. 
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